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^ ' Cluster plasmas are magnetised already at very low magnetic field strength. Low collision- 

ality implies that conservation of the first adiabatic invariant results in an anisotropic viscous 
stress (Braginskii viscosity) or, equivalently, anisotropic plasma pressure. This triggers fire- 
hose and mirror instabilities, which have growth rates proportional to the wavenumber down 
to scales of the order of ion Larmor radius. This means that MHD equations with Braginskii 
viscosity are not well posed and fully kinetic description is necessary. In this paper, we review 
the basic picture of small-scale dynamo in the cluster plasma and attempt to reconcile it with 
the existence of plasma instabilities at collisionless scales. 
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Abstract 



The intracluster medium (ICM) is a unique environment for studying astrophysical tur- 
bulence: it is made of amorphous diffuse fully ionised plasma with almost no net rotation or 
shear. The ICM is, therefore, in a certain sense, a case of "pure" turbulence. The presence 
t— I | of magnetic fields in the ICM with B ~ 1...10 /iG and field scale ~ 1 kpc is well documented 
(e.g., |Carilli fc Taylor 2002] ). Because of the availability of the rotation measure data from ex- 
tended radio sources, it is possible to infer the magnetic-energy spectra for some of the clusters 



■ with spatial resolution of - 0.1 kpc jVogt fc Enfilin 2003[|Vogt et al. 2004|). This is due t o 
Or improve dramatically with the arrival of LOFAR and, especially. SKA l|Gaensler et al. 20f)4J) . 

At the same time, there is a rising interest among astronomers in measuring the velocity spec- 
?H ' tra of the cluster turbulence using instruments both future (ASTRO-E2, |Sunyaev et al. 2 003) 
and present (XMM-Newton, Sch uecker et al. 2004J) . Once we have both kinetic and magnetic- 
energy spectra, a clear observational picture of the magnetised cluster turbulence will emerge. 
The purpose of this paper is to review our theoretical understanding of this turbulence. 

There is no consensus on the exact mechanism that drives turbulence in clusters. The 
main candidates are galaxy wakes, jets, and cluster mergers. We will leave a discussion of 
these outside the scope of this paper. The important point for us is that all these mechanisms 
imply the outer scale of the cluster turbulence in the range Iq ~ 10 2 ...10 3 kpc with outer scale 
velocities comparable to the thermal (sound) speed vth = (T/nii) 1 ^ 2 ~ 10 3 km/s, where m, 
is ion mass (assuming hydrogen plasma) and T ~ 10 8 K the plasma temperature (assuming 
for simplicity Tj = T e ). This gives the outer eddy timescale of order l^/vth ~ 10 9 yr. 

The ion kinematic viscosity is v ~ viu/va, where va = Airne 4 In Am i ^j 1-3 / 2 i s the ion- 
ion collision frequency, n ~ 10 _2 ...10~ 3 cm -3 is the ion number density (for cluster cores), 
e is the electron charge, and In A ~ 20 is the Coulomb logarithm (Spitzer 1962]). This gives 
Reynolds numbers in the range Re ~ 10 2 ...10 3 . As the turbulent velocities at the outer scales 
are ~ Vth, the turbulence in the inertial range will be subsonic. It is natural to assume that 
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Kolmogorov's dimensional theory should apply at least approximately. Then the viscous scale 
of the turbulence is l v ~ loRe~ 3 ^ 4 ~ 10. ..30 kpc. These numbers appear to agree quite well 
with observations of turbulence in the Coma cluster via pressure maps ( Schue ck"er et al. 20040 . 

Because the ICM plasma is fully ionised, the viscous cutoff is, in fact, quite different from 
the usual Laplacian viscosity that normally appears in MHD equations. It can be shown 
that at frequencies below the ion cyclotron frequency f2^ = eB/crrii and scales above the ion 
Larmor radius pi — v t \JSli, fluid velocity and magnetic field satisfy 



du B 2 \ „ r~, N i B-VB 

P Tt = -V^ + ^J+V.[bb(px-p„)] + - ir - l (1) 

— = B ■ Vu- BV • u + ?7V 2 B, (2) 

where d/dt — d/dt + u • V, p = rmn is mass density, b = H/B, p± and pu are the plasma 

pressures perpendicular and parallel to the magnetic field, r\ ~ T _3 / 2 m;|/ 2 e 2 c 2 In A/ An is the 
magnetic diffusivity (S pitzer 1962| we ignore the order-unity difference between rj^ and r/±). 
Eq. Q is the appropriate description for the fluid motions in the ICM provided the ions are 
magnetised, i.e., pi <C A m f p , where A m f p ~ Vth/^ij ~ 1...10 kpc is the ion mean free path. 
This requirement is satisfied if B ^S> 10~ 18 G, which is far below the observed field strengths 
of 1...10 pG. Note that 10 pG is roughly the field strength corresponding to the energy of 
the viscous-scale eddies (~ pf 2 h Re -1 ^ 2 ), which is the lower bound for dynamically important 
fields. Thus, even for dynamically weak fields, the plasma is already very well magnetised. 

The plasma pressures in Eq. Q must be calculated kinetically. At low frequencies (u> <C 
kvth) and at scales above A m f p , the total isotropic pressure p± + B 2 /8tt can be determined 
from incompressibility, V ■ u = 0, while the pressure anisotropy is calculated perturbatively 
from the kinetic theory ( |Braginskii 1965| ): 

pj_ - p|| = ipv bb : Vu ee pv 2 h A, (3) 

where v ~ WthA m f p ~ v 2 h /vu is the ion viscosity and A is the dimensionless measure of 
the pressure anisotropy. The pressure anisotropy is proportional to (the parallel component 
of) the rate of strain, so it is determined by the viscous-scale eddies' turnover rate u v /l u . 
For Kolmogorov turbulence with outer-scale velocity ~ Dth, we have u„/i>th ~ Re -1 / 4 and 
k/Amfp ~ Re~ 1/4 , so A - UvXjofp/vfalv ~ Re _1/2 . 

Eqs. (|1I3|I plus the incompressibility condition V • u = are a closed system that we 
will call the Braginskii MHD. Two key properties of these equations should be noted. First, 
velocity gradients transverse to the magnetic field (e.g., shear- Alfven-polarised fluctuations) 
are not dissipated. Velocity fluctuations can now penetrate below the viscous cutoff. Second, 
the velocities that are dissipated are those that change the strength of the magnetic field: 
indeed, from Eq. Q), 

1 dB C£ 

B^- hh - VU - (4) 

Physically, the Braginskii viscosity encodes a fundamental property of charged particles mov- 
ing in a magnetic field: the conservation of the first adiabatic invariant p = mv 2 _/2B, which is 
only weakly broken when A m f p 3> pi. This means that any change in B must be accompanied 
by a proportional change in p±. Therefore, the emergence of the pressure anisotropy is a 
natural consequence of the changes in the magnetic-field strength and vice versa. 
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Let us ignore for a moment the inability of the Braginskii viscosity to damp all velocity 
fluctuations and consider the turbulent motions between the outer scale and the viscous 
scale and the random stretching of the magnetic field by these motions, or the small-scale 
dynamo. Such a view might appear justified because the Braginskii viscosity cuts off any 
motions that change the field strength, so the small-scale dynamo action is associated only 
with velocities above the viscous scale. The turbulent small-scale dynamo is a problem that 
we have studied extensively (for the full account of our work and a list of relevant references, 
see Schckochihin et al. 2004). Here we reiterate briefly the main points. 

In the weak-field regime (no back reaction on the flow), the turbulent stretching is done 
primarily by the viscous-scale eddies because they turn over the fastest. The magnetic energy 
grows exponentially roughly at their turnover rate u v /l v ~ 10~ 8 yr _1 . The turbulence tends 
to arrange the magnetic fields in long thin flux sheets (or ribbons): the field reverses its 
direction at the smallest scale Ib available to it (if the magnetic cutoff is modelled by the 
Ohmic diffusivity n, it is the resistive scale: Ib ~ l v ~ l v VxoT x ^ , where Pm = v/rj) but field 
lines curve at the scale of the flow (L ~ l v , the viscous scale) except in the sharp bends of 
the folds. The field is in the state of reduced tension everywhere: both in straight and bent 
parts of the folds, BK 1 / 2 ~ const, where K = |b • Vb| is the field-line curvature. Thus, the 
field strength and the field-line curvature are anticorrelated. 

The tendency to fold the field is a kinematic property of random stretching. A key 
numerical finding is that the folded structure is also a feature of the nonlinear saturated state 
of the dynamo. The folded structure allows the small-scale direction-reversing magnetic field 
to back react on the flow in a spatially coherent way. The back-reaction mechanism we have 
proposed consists in the velocity gradients becoming locally anisotropic with respect to the 
direction of the folds (with bb : Vu partially suppressed) so that the dynamo saturates at the 
marginally stable balance between reduced "parallel" stretching and "perpendicular" mixing. 

In the kinematic regime, the characteristic fold length is the viscous scale, lu ~ l v . In the 
saturated state, the folds tend to elongate to the outer scale, In ~ ?o- We have proposed that 
the saturation is controlled by the anisotropisation of the outer-scale eddies by the folded 
fields, while the inertial range (Iq > I > l u ) is populated by Alfven waves that propagate 
along the folds with the dispersion relation uj 2 = kk : hh(B 2 )/4:irp. 

The physical considerations outlined above provide a qualitative framework for interpret- 
ing the saturated spectra of magnetic and kinetic energies that are seen in the numerical 
simulations of forced isotropic MHD turbulence: magnetic energy significantly exceeds ki- 
netic energy at small scales and, in the case of large Pm, tends to concentrate below the 
viscous scale (see, e.g,. ISchekochihin et al. 20041 |Haugen et al. 2 004). Such a magnetic spec- 
trum makes sense if it is thought of as a superposition of folds with field reversals at the 
resistive scale and and Alfven waves propagating along the folds. 

Recently published magnetic-energy spectra based on RM measurements from extended 
radio sources fo r Abell 400, 2255, 2634, and Hydra A clusters (see |Vogt fc Enfilin 2003| 
|Vogt et al. 2004[ and Enfilin's contribution in these Proceedings) show magnetic energy to be 
peaked at ~ 1 kpc with the overall shape of the spectra quite similar to those that emerge in 
numerical simulations cited above. Does this mean that these simulations correctly describe 
cluster turbulence? The answer to this question is probably in the negative. All existing sim- 
ulations cither use isotropic viscosity and Ohmic diffusivity to model kinetic and magnetic 
cutoffs or rely on numerical dissipation. Clearly, simulations with a magnetic cutoff at or not 
far below 1 kpc (i.e., about a decade below the viscous scale) are likely to produce magnetic 
spectra that will look roughly like the observed ones. However, if we estimate the resistive 
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cutoff scale based on the formula for Ohmic diffusivity given above, we will get l v ~ 10 4 km 
— many orders of magnitude below 1 kpc. Furthermore, no published numerical study has 
included the effect of anisotropic viscosity to address the undamped velocity fluctuations at 
subviscous scales. In what follows, we show that the Braginskii viscosity not only fails to 
efficiently damp kinetic energy at the viscous scale but triggers fast-growing instabilities at 
subviscous scales. This fundamentally alters the structure of the turbulence at these scales. 

It has, in fact, been known for a long time that pressure anisotropy leads to plasma instabil- 
ities IjRosenbluth 19561 lUhandrasekhar et al. 19581 IParker 19581 [VeH enov fc Sagdeev 1958| ). 
A vast literature exists on these instabilities, which we do not attempt to review. The stan- 
dard approach is to postulate a bi-Maxwcllian equilibrium distribution with Tj_ 7^ Tji (see, 
e.g.. lFerriere fc Andre 20Q 2land references therein). We do not need to adopt such a descrip- 
tion because Eqs. (Ill- ill incorporate the pressure anisotropy in a self-consistent way. 

Consider a turbulent cascade that originates from the large-scale driving, extends down 
to the viscous scale, and gives rise to velocity and magnetic fields u and B that vary on time 
scales ~ (Vu) _1 and on spatial scales ~ l v . We study the stability of such fields. The presence 
of turbulent shear (velocity gradients) gives rise to the pressure anisotropy given by Eq. J3J • 
We now look for linear perturbations <5u, <5B, Spx, Sp« that have frequencies u> 3> Vw, and 
wavenumbers k ^> l~ x . With respect to these perturbations, the unperturbed rate-of-strain 
tensor Vu can be viewed as constant in space and time. Linearising Eq. Q and neglecting 
temporal and spatial derivatives of the unperturbed quantities, we get 



itupSu = -ikS ( px + ^—j + (px -p\\ + ^—\ S(b ■ Vb) 



b<5 



8ir J V 11 47r 

B 2 \ V||B 



V|| (px -P\\) - [px-P 



4vr J B 



(5) 



and, from Eq. J2J), Sh = — (fcii/w)<Ju_L and 5{h- Vb) = iknSh. We have denoted V|| = b- V and 
fen = b • k. In the resulting dispersion relation, it is always possible to split off the part that 
corresponds to the modes that have shear-Alfven-wave polarisation, 8u\\ = 0, kj^ • Su± = 0: 

c 2 = fc[|i4(A + 2/r 1 ), (6) 

where (3 = 2v^/v\ = 8iTpv 2 h / B 2 . When the magnetic energy is larger than the energy of 
the viscous-scale eddies, (3 <C |A| _1 ~ Re 1 ^ 2 , Eq. © describes shear Alfven waves. In the 
opposite limit, (3 3> |A| _1 , an instability appears: when A < 0, the growth rate is 7 = 
fcllUthlAI 1 / 2 ~ fc||W t hR.e -1 ^ 4 . Since fcii 3> the instability is faster than the rate of strain, 
7 3> Vw ~ WthRe -1 ^ 4 //^, in accordance with the assumption made in the derivation of Eq. 10. 
This instability, triggered by p\\ > px, is called the firehose instability. It does not depend 
on the way the pressure perturbations are calculated because it arises from the perturbation 
of the field-line curvature b • Vb in Eq. J5J : to linear order, it entails no perturbation of the 
field strength [from Eq. J5J, SB/B = — (k\\/u})8u\\\ and, therefore, does not alter the pressure. 

In order to determine the stability of perturbations with other polarisations, 5px and 
5p\\ have to be computed. We consider scales smaller than the mean free path, fcA m f p ^> 1, 
where the plasma is collisionless. It is sufficient to look for "subsonic" perturbations such that 
lu <^ kvth because the high-frequency perturbations are subject to strong collisionless damping 
and cannot be rendered unstable by a small anisotropy. A straightforward linear kinetic 
calculation ( Schekochi hin et al. 2005|) shows that the slow-wave-polarised perturbations are 
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unstable with the growth rate 
/2\ 1/2 

Modes with fen > \/2kj_ are unstable if A < (slow- wave-polarised firehose instability); modes 
with k± > kii/y/2 are unstable if A > (mirror instability). 

The instabilities we have described have growth rates oc fen. Thus, MHD equations with 
Braginskii viscosity do not constitute a well-posed problem. This means, for example, that 
any numerical simulation of these equations will blow up at the grid scale unless additional 
isotropic viscosity is introduced (this has been confirmed by J. L. Maron, unpublished). 

Eqs. © and Q are valid in the limit kpi — ^ 0. If finite pi is taken into account, the 
growth rates peak at fc ~ pj l with 7 max ~ |A|f2j for the mirror and 7 max ~ lAI 1 / 2 ^ 
for the shear-Alfven-polarised firehose. For clusters, this gives values of 7 max in the range 
1CP 8 ...10~ 7 (-B/1CP 18 G) yr _1 — much faster than the viscous-eddy turnover rate. 

We have shown that, given sufficiently high 0, the firehose and mirror instabilities occur in 
the regions of decreasing (A < 0) and increasing (A > 0) magnetic-field strength, respectively. 
The small-scale dynamo action by the viscous-scale motions associated with the turbulent 
cascade from large scales will produce regions of both types: the folded structure explained 
above contains growing straight direction-reversing fields and weakening curved fields in the 
corners of the folds. This structure is intrinsically unstable: straight growing fields to the 
mirror, curved weakening fields to the firehose instability. It seems plausible, however, that 
the magnetic field does grow until /3 ~ | A| _1 ~ Re -1 / 2 , so the instabilities are stabilised. This 
corresponds to B ~ 10 p,G, which is comparable to the observed field strength in clusters. 

After the instabilities are quenched, the fluid motions can be separated into two classes: 
above the viscous scale, there is the usual turbulent cascade, to which the picture of small- 
scale dynamo given in Schckochihin et al. 2004 and outlined above applies; below the viscous 
scale, the slow-wave-polarised motions are damped, while the shear Alfven waves (the stable 
counterpart of the firehose instability) can exist. If we conjecture the survival of the folded 
structure, which is a generic outcome of random stretching, then the Alfven waves will prop- 
agate along the folds. Note that the linear stability analysis outlined above can be recast in 
terms of the perturbations of striped direction-reversing fields represented by the tensor BB. 
What we still lack is an estimate of the reversal scale. 

Even for rms field strengths sufficient to quench the instabilities, there will always be 
regions where the field is weak. These regions can become unstable unless they are smaller 
than the local ion Larmor radius. We conjecture that the overall structure of the magnetic field 
is determined by the requirement that the field scale in the weak-field regions should be smaller 
than local value of pi. The following rather tentative argument is an attempt to estimate 
the characteristic field scale based on this idea. As we mentioned above, the folded structure 
is characterised by the anticorrelation between the field strength and field-line curvature: 
BK 1 / 2 ~ const. The weak-field regions are the bending regions where the field is curved (see 
Fig. [Q. A simple flux-conservation argument gives B\, erL ^/B Trns ~ Ib/1\\, where B rms is the 
field in the strong-field regions, Ib is the reversal scale, and Zy is the fold length, which, in 
the nonlinear regime, can be as long as the outer scale Iq of the turbulence. We estimate 
the curvature in the bending regions by ifbend ~ 1/Pi,bcnd, where p^bend = PiB lms / B hcnd is 
the ion Larmor radius in the bending regions and p; = i>th/^i = cm^T 1 / 2 / 'eB rms is the ion 
Larmor radius in the strong- field regions. The curvature in the strong-field regions is ~ . 
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Figure 1: A sketch of the folded structure and scales referred to in the derivation of Eq. ©. 

Then -Bbcnd/-E>rms ~ (/Oi.bcnd/'n) 1 ^ 2 - Assembling these relations and substituting numbers, we 
obtain the following estimate 

This is still substantially smaller than the scale at which the observed spectrum peaks, but 
not entirely unreasonable, considering the simplistic nature of our argument. Note that the 
shear Alfven waves can exist at scales below Ib, with a cutoff at k±pi ~ 1, so the turbulence 
spectra are likely to have extended power tails. 

A somewhat more detailed exposition can be found in ISchekochih in et al. 2 005 
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